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Abstract. Adapting a method developed for the study of quantum chaos on quantum 
(metric) graphs [1], spectral C functions and trace formulae for discrete Laplacians on 
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CO ' class of C. functions proposed originally by Ihara [2\, and expanded by subsequent 
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authors [31 3]. Finally, a model of "classical dynamics" on the discrete graph is 
proposed. It is analogous to the corresponding classical dynamics derived for quantum 
graphs [I]. 



1. Introduction and preliminaries 

Some ten years ago quantum (metric) graphs were proposed as a convenient paradigm 
for the study of quantum chaos in compact [T] and scattering [5] systems. The 
crucial point which highlighted the close similarity between metric graphs - whose only 
claim to complexity is their topology - and chaotic Hamiltonian flows is the formal 
similarity between the trace formulae O [7] which express the spectral densities as 
sums over periodic orbits. Requiring additionally that the lengths of the bonds are 
rationally independent, and that the graph is well connected, render the spectrum of 
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the Schrodinger operator on graphs sufficiently disordered to display spectral statistics 
which are consistent with the predictions of Random Matrix Theory. Another important 
contact point was the identification of the classical dynamics on the graph, which is 
derived from the quantum evolution, and which describes, under certain conditions, an 
exponential approach to equilibrium, in analogy with mixing Hamiltonian fiows. The 
( function for the Peron-Frobenius operator on the graph, can be written in terms 
of the same periodic orbits which are used for the quantum spectral ( function - in 
close similarity with the expansion of the Ruelle ( function as a sum over periodic 
orbits for chaotic Hamiltonian fiows. The derivation of the trace formula for quantum 
graphs which was presented in |lj differs from the original method [6]. It uses another 
approach (called sometimes the "scattering approach"), which reveals in a natural way 
the underlying classical dynamics. 

Discrete graphs, where only the graph topology and not its metric plays a role 
are mostly studied in number theory, combinatorics etc. There is abundant literature 
relating to various aspects of graphs. Much of the relevant material to the present 
discussion can be found in [H |9]. Audrey Terras' review [10] surveys the field, and 
its relation to quantum chaos. The present work attempts to highlight further this 
quantum chaos connection, by proposing trace formulae and spectral ( functions, and 
linking them with the Ihara ( function [2] and some of its recent generalizations. To 
introduce these concepts, a few preliminaries and definitions are necessary, and they are 
provided below. 

A graph Q consists of V vertices connected by B bonds. The V x V connectivity 
(or adjacency) matrix C is defined such that Ci,j = 1(0) if the vertices i,j are connected 
(disconnected). Graphs with parallel bonds or loops are excluded. The valency (some 
times referred to as the degree) of a vertex is the number of bonds which emanate from 
a vertex. It is denoted by Vi = YlJ=iCi,j- To any bond b = {i,j) one can assign an 
arbitrary direction, resulting in two directed bonds, d = {i,j) and d = {j,i)- Thus, 
the graph can be viewed as V vertices connected by bonds b = 1, ■ ■ ■ , i? or by 2B 
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directed bonds d = 1, - ■ ■ , 2B. (The notation b for bonds and d for directed bonds will 
be kept throughout). It is convenient to associate with each directed bond d = 
its origin o{d) = i and terminus t{d) = j so that d points from the vertex i to the 
vertex j. The bond d' follows d if t{d) = o{d'). A periodic orbit (cycle) of length n is a 
sequence of n successively following directed bonds di, - ■ ■ ,dn and di follows (i„. Cyclic 
permutations of the bonds generate the same periodic orbit. A primitive periodic orbit 
is an orbit which cannot be written as a repetition of a shorter periodic orbit. The set of 
primitive n-periodic orbits will be denoted by V{n), and V = IJ^2 ^('^)- An important 
subset of V{n) is the set of n primitive periodic orbits without back-scatter, namely, 
periodic orbits where di^i ^ di. The corresponding sets will be denoted by C(n) and 

c = u:r=2C(^)- 

The Laplacian of a discrete graph is defined as 

L = -C^D, (1) 

where C is the connectivity matrix, and D is a diagonal matrix with Di i = fj. It is 
a self-adjoint operator whose spectrum consists of V non negative real numbers. The 
spectrum is determined as the zeros of the secular function (characteristic polynomial) 

Zi(A) = det(A/(^)-L) . (2) 

Here, A is the spectral parameter and /*^^^ is the unit matrix in V dimensions. The 
lowest eigenvalue is 0, and it is simple if and only if the graph is connected. 

It is sometimes convenient to generalize the Laplacian ([1]) by replacing the matrix 
C by a matrix C whose zero entries coincide with those of C, but arbitrary, strictly 
positive weights Wjj- (= w^^i) replace the values 1 where Cjj = 1. One then defines 
Di i = Ui = Ci j and the generalized Laplacian is 

L = -C + iD. (3) 

The spectrum of L consists of the zeros of the secular equation (characteristic 
polynomial) Zi{\) = det(A/*^^^ — L). The spectrum is non negative, is in the spectrum 
and it is a simple eigenvalue if and only if the graph is connected. 
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The focus of the present work is on ( functions and trace formulae for discrete 
graphs. This research subject was initiated by Ihara [2j who defined a ( function for a 
graph as 

c(«r =n(i-«")"^"^ (4) 

n 

where |C(?t,)| is the cardinahty of the set C{n), and m G C with |m| sufficiently small to 
ensure the convergence of the infinite product. Following Ihara's original work, several 
authors (see e.g., [Tl] for a survey of the methods) have proved that 

C{u)-^ = (1 - u^y-^ det(j(^) -uC + u^Q) . (5) 

Here, r = B — V + 1 is the rank of the graph (the number of independent cycles on 
the graph or equivalently, the rank of its fundamental group). J*-^-' is the unit matrix 
in V dimensions, C is the connectivity matrix, and the diagonal matrix Q = D — 
If the graph is w-regular, that is Vi = v Vi, the non trivial poles of the Ihara ( (the 
trivial poles are at u = ±1) can be easily computed from the eigenvalues of the graph 
Laplacian ([T]). 

The following ( function defined by H.M. Stark [12] will serve as an example of the 
more recent developments in the field. Consider a matrix Y in the space of directed 
bonds 

Yd',d = Vd',d 5o(d'),tid) (1 - 5rf/,rf)- (6) 

where rjd'^d are arbitrary. Note that matrix elements between reversed bonds are 
excluded. Associate with any primitive periodic orbit c G C the amplitude 

fc = Vdn,dn-1 Vd„-l,dn-21 ' ' ' Vd2,dl Vdi,dn ■ (7) 

Then, 

CEiY)-' ^ J] (1 - /,) = det(/(2^) - Y), (8) 

where /(^^^ is the unit matrix in 2B dimensions. This result will be used in the last 
section. 
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In the next section, other ( functions are defined, discussed, and expressed as 
rational functions which are reminiscent of ([5]) and ([8]), but are different in many 
respects. Trace formulae for the spectra of the Laplacians (fT][3l) will also be derived. 
In the last section, the approach developed here will be compared with its analogues in 
the theory of quantum graphs, and the "classical dynamics" on the discrete graph will 
be proposed. 

2. Secular functions, ( functions and trace formulae 

To start, an alternative form of the secular equations for the Laplacians (IDISll will be 
derived. It is convenient to begin with a detailed derivation for the traditional Laplacian 
([1]). The necessary modifications for the generalized form will be indicated later. For 
both Laplacians, the secular function will be shown to take the form 



where U{X) is a unitary matrix of dimension 2B which depends on the spectral parameter 
A. By construction, Zs{X) is real for A G M, and its zeros will be shown to coincide 
(with their multiplicity) with the spectrum of the Laplacian. Thus Zs{\) and Zl{\) can 
differ at most by a multiplicative function of A which does not vanish for real A. This 
construction of the secular function paraphrases the "scattering approach" introduced 
in [T] for quantum graphs. (Derivations which are similar in spirit were discussed in 
[T3t [H] , see also [I5] and references cited therein) . 

To compute an eigenvector ip = {ipi, ■ ■ ■ ,ipv) of L, corresponding to an eigenvalue 
A, the following steps are taken. To each bond b = one associates a bond wave 
function 



Zs{\) = ^ (det t/(A))-^ det (/(^^) - f/(A)) 



(9) 



X G {±1} 



(10) 



subject to the condition 



(11) 
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Consider any vertex indexed by i, and the bonds (61, 62, •••^i>J which emanate from i. 
The corresponding bond wave functions have to satisfy three requirements in order to 
form a proper eigenvector of L. 

I. Uniqueness: The value of the eigenvector at the vertex i, ipi, computed in terms of 
the bond wave functions is the same for all the bonds emanating from i. The following 
Vi — 1 independent equalities express this requirement. 

aft, e*? + ag^ e"*^ = a^^ e*^ + ag^ e"*? = . . . = a^^. e*^ + e"*^ . (12) 

//. is an eigenvector of L : At the vertex i, X^jLi = terms of the bond 

wave functions this reads, 

- J2 [«^. e-^f + a^, e+^f ] ={X-v.)l^J2 i^"- + e'^^] . (13) 

1=1 ' m=l 

To get the equation above, ipi was presented as 

Together, fll2p and fll3p provide fj homogeneous linear relations between the 2vi 
coefficients a^, where d stand for directed bonds which are either incoming to {t{d) = i) 
or outgoing from {o{d) = i) the vertex i. Using these equations, the outgoing coefficients 
are expressed in terms of the incoming ones. 



ad= ^ (^%{^) o-d' Vd : o{d) = i , (15) 

d' : t{d')=i 



where. 



_,iai(A) 



l + i{l- ^) 



(16) 



The vertex scattering matrices a^^\\) are the main building blocks of the present 
approach. They distinguish clearly between back-scatter transitions {d = d') and the 
transitions to other bonds, for which the same strength is given, independently of the 
original and the final bonds. For real A the vertex scattering matrices are unitary 
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matrices and they are the discrete analogues of the vertex scattering matrices derived 
for the Schrodinger equation on graphs [1]. 

///. Consistency : The hnear relation between the incoming and the outgoing coefficients 
( TT5l) must be satisfied simultaneously at all the vertices. However, a directed bond {i,j) 
when observed from the vertex j is outgoing, while when observed from i it is incoming. 
This consistency requirement is implemented by introducing the Evolution Operator 
Ud'.dW in the 2B dimensional space of directed bonds, 

t/.',.(A) = 5*(,),oK)^2!f(A). (17) 

{U is also referred to in the literature as the Bond Scattering Matrix [Ij). The evolution 
operator is unitary U U'^ = /(^-^^ for A G M due to the unitarity of its constituents cr^*\ 
Denoting by a the 2B dimensional vector of the directed bonds coefficients defined 
above, the consistency requirement reduces to, 

f/(A)a = a. (18) 

This can be satisfied only for those values of A for which 

e(A) = det (J(2^) - U{\)) = . (19) 

For real A the spectrum of U{X) is restricted to the unit circle. Therefore |^(A)| is finite 
for all A G R. Due to (fT6|) the matrix elements of U{\) are ratios of monomials in A. 
These two properties imply that ^(A) = p(A)/g(A) where p and q are polynomials of 
the same degree in A, and their degree is at most 2B. The zeros of g(A) coincide with 
the poles of det U{\). They are complex because | det ?7(A)| = 1 for A G M. A straight 
forward computation yields, 

det^(A) = n i_,n -I) ' ^ q{X) = Const W ( 1-^(1--)) • (20) 
j=i H j=i ^ ^ 

Thus, det U has exactly V complex poles, implying that the degree of p(A) which equals 
the degree of g(A) is also V. Note finally that the zeros of p{X) coincide with the zeros 
of the secular function Zl(\) = det(A/*^^^ — L) which is also a polynomial of degree 
V. Hence, p{X) and Zl{X) are identical up to a constant factor. It is convenient to 
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define the secular equation so that it is real on the real axis. This can be achieved 
by multiplying ^(A) by (det f/(A))~5. A further factor of normalizes the resulting 
function to approach 1 as |A| — oo. The resulting secular equation reads 

Zs{X) = ^ (det f/(A))-^ det (J^^^) - U{X)) (21) 

J_yr A+^(l-t) y p(A) ^ det (A/(^) - L) 

This expression for the secular equation is the basis for the further results of the present 
work. To begin, use is made of the fact that the spectrum of U{X) for Xm(A) < is 
confined to the interior of the unit circle. Thus, for for any A with an arbitrarily small 
(but finite) negative imaginary part, we expand 

oo _ 

logdet(/(2^) - U(X)) = - y - trt/"(A) , (22) 

and 



n 

n=l 



trf/"(A) = J2 ^ Yl "P^^) • (23) 

m:m\n pG'P{m) 

The sum above is over all the primitive periodic orbits p with period m which is a divisor 
of n, p = di, ■ ■ ■ ,dm and 

ap{X) = (Td,4^ (A) ■ ■ ■ ad.,,d, (A) . (24) 

The explicit dependence of ap(A) on A is obtained from the following expressions for the 
vertex scattering matrix elements, 

e for d' ^ d, 

2{v.-\) (25) 



4(fi-l) 



1 



1- 



1 

2 —I arctan r,— r^- 



e (.,-i)^+(.,~A)^-i fo^ d' = d 



where j = t{d) = o{d'). The explicit expressions above were written so that for real A 
the absolute square of the is a product of "transition probabilities" , while the phase of 
Qp is a sum which plays the role of the "action" or "length" associated with the periodic 
orbit. Substituting fl23|) in fl22l) . and summing over the repetition numbers ^ one gets, 

det(j(2^) - U{X)) = - ap{X)) . (26) 
pev 
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The ( function which is introduced in the present work is defined as 

Cs{X)-' = ll(^-a,{X)). (27) 

pev 

Combining ( 122|) and ( l26l) with the definition of Cs{X) gives 

n;.i + - A)) 

This is one of the main results of the present work. It provides a "Ihara" - hke identity 
which expresses an infinite product over primitive periodic orbits on the graph in terms 
of the characteristic polynomial of the graph discrete Laplacian. The main difference 
is that here, all the periodic orbits, including orbits with back-scattering, contribute to 
the product, and that the amplitudes Op depend on the spectral parameter in a more 
complicated way. To get a closer look at the (s function and its relation to the Ihara ( 
function, it is instructive to write (s for a general f-regular graph. For this purpose, it 
is convenient to define a new complex variable, 

l+i(l--) , , 

^ = IV 29 

which is a 1 1 map of M to the unit circle in C. With these simplifications, Cs{^) fl28l) 



reduces to 



2z \ , / _ .2 — 1 



It is convenient to define 'ysiz) = z^(siz), in terms of which a functional equation for 
(s can be written as 

^s{z-') = ilsizir , (31) 

where (■)* stands for complex conjugation. Functional equations of similar type 
are satisfied also by the Ihara ( function (for f-regular graphs) as well as by most 
other functions of this genre. Typically, functional equations enable the analytical 
continuation of ( functions which are defined by infinite products, beyond their radius 
of convergence. Here also it provides the analytic continuation of Cs{^) to the exterior 
of the unit disc. 
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The periodic amplitudes ap{z) simplify considerably for f-regular graphs. Denote 
by Up the period of the primitive periodic orbit p, and by jSp the number of vertices in 
p where back-scattering occurs: /?p = jj |i : di = di+i, di G p, dn^+i = c^ij- Then, 

apiz) = e-i"^ ( ^ V' (-1)^^ fl - (32) 



The results above pave the way to the derivation of trace formulae for the discrete 
Laplacians. Trace formulae provide a powerful tool in spectral theory. They express 
the spectral density (written down formally as a sum of Dirac 6 functions located at the 
spectral set) in terms of information derived from the manifold metric. The spectral 
density is written as a sum of two contributions - both of which have a geometric origin. 
The first is a smooth function of A whose asymptotic limit at A — oo was first studied 
by Weyl. The second contribution is an infinite sum over periodic geodesies on the 
manifold. The equality between the spectral density and its geometric representation 
should be understood only in the sense of distributions. An analogous trace formula will 
be derived now for the discrete Laplacian. Making use of Cauchy theorem and the fact 
that Zs{X) is analytic in the vicinity of the real A axis, and real on it, one can write, 

V 

d{X) = ^5(A- Aj) (33) 
= - lim Xm^r log Zs(X - ie) . (34) 

71 e^0+ dA 

Using 

Zs{\) = ^(detf/(A))-^det(/(2^)-t/(A)) , (35) 
the explicit form of detf/(A) ( 120|) and the periodic orbit expansion ( l23l) . one gets, 

j=l J ^ Vj' r=l pg-p ^-^^ 

The first term is the "smooth" (Weyl) contribution to the spectral density. It consists 
of a sum of Lorenzians with poles at Xj = Vj{l ± i). This sum is analogous to Wigner's 
semi-circle density in Random Matrix Theory. 
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The explicit expression for the fluctuating part can be written down explicitly using 

( I25l) . Noting that the ap(A) are complex numbers with A dependent phases 0p(A), the 

periodic orbit sum in the trace formula is a fluctuating function of A. It is the term 

which turns the r.h.s. of ( l36l) to a distribution when e ^ 0. 

So far, the discussion was restricted to the "traditional" Laplacians. The extension 

to the generalized Laplacians, starts by modifying the definition of the bond wave 

functions (ITOl) to read, 

^6 = v^6(afee*^" + a^e-*^^) . (37) 

Then, following the same steps as above, the vertex scattering matrices are derived, and 
they take the form 

^&(A)=^(5.>-^(l + e^"'(^Mv/^) ; e--W= j^^ll^fj , (38) 

where Uj = Wij as defined previously. The subsequent derivation follows the same 
steps, resulting in the generalized (s function, 

C,(A)- . Ui'-a^iX)) = . (39) 

A trace formula is also derived in the same way, 

^ ~ 1 d 

d{X) = ^5(A- Aj) = - lim Jm— logZ5(A-ie) (40) 

i=i 

= -y- — ^—T^ - -^^^y y -i^M e^^^^^^^ 

7T ^ Ui 1 + (1 - TT dA ^ ^ n(p) ' ^' 

j=l ■' ^ Uj ' r=l peV ^ ' 

The expressions for ap{X) can be derived by a simple modification of ( |25l) and therefore 
they will not be written down here. 

3. Classical dynamics 

The present approach emerges from the alternative secular function for the spectrum 
of Laplacians, based on the quantum evolution operator U (A) in the space of directed 
bond amplitudes a G /^(C^"^). Consider U{X) as a quantum map which maps this 2B 
dimensional space onto itself. U is unitary and hence the map conserves the P norm - 
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the quantum probability. The condition U{X)a = a can be interpreted as a requirement 
that A„ is an eigenvalue if there exists a non trivial vector a which is stationary under the 
action of the quantum map [16j . The requirement of stationarity is naturally associated 
with the eigenvalue being in the spectrum of the underlying Hamiltonian. 

The building blocks for the theory are the vertex scattering matrices. Similar 
matrices appear in the theory of quantum graphs. There, they emerge when the 
Schrodinger equation on the graph is augmented by vertex boundary conditions which 
render the resulting operator self adjoint. The self adjoint extension is not unique, and 
depends the spectral parameter k and on an arbitrary parameter k, which interpolates 
between the "Dirichlet" (k = 0) and the "Neumann" {k = oo) boundary conditions 
[U |T7] . The scattering matrices for discrete graphs are obtained from their quantum 
graph analogues by replacing n/k by A. 

The unitary quantum evolution operator is the starting point for the construction 
of a classical evolution on the discrete graph. The classical "phase space" in this case 
are the probability vectors p G P^R."^^) where the components of p are interpreted 
as the probabilities to find the classical system on the corresponding directed bonds. 
The classical transition matrix is constructed from the quantum probability to make a 
transition from d to d' 



The unitarity of U implies that M is bi-stochastic, namely, ^^^^'4 = ^d',d = 1- 
This transition matrix induces a discrete, random walk dynamics in phase space. If n 
denotes the discrete "time" , 



This Markovian evolution preserves the norm - the classical probability. The spectrum 
of M is confined to the interior of the unit circle. 1 is always an eigenvalue corresponding 
to an eigenvector with equal components which describes the system in an equilibrated 
state. When the eigenvalue 1 is the only eigenvalue on the unit circle, the classical 



Md'd 




(41) 



p(n + 1) = Mp{n) . 



(42) 
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dynamics drives the system to equilibrium at a rate which depends on the distance of 
the next highest eigenvalue to the unit circle. This classical dynamics is identical to the 
one which was introduced in the study of quantum graphs [Ij. It plays an important 
role in the theory of spectral statistics on quantum graphs [H [TBI [191 EO]. Finally, the 
analogue of the Ruelle ( function for the evolution induced by P can be easily written 
down starting with the secular function 

Zm(/x) = det(/(2^) - /iM) (43) 

The periodic orbit sum is identical to the trace formula (136|) in which the amplitudes 
are replaced by their absolute squares. 

To emphasize the intricate connections between the concepts developed here and 
their predecessors |2ll3lll], consider a f -regular graph (t> > 2), and the classical evolution 
operator obtain for the spectral parameter A = v + i{v — 2) , corresponding to z = v — 1 
in (129|) . At this value, o"^^' = and ad' 4 = 1 for d' 7^ d. The resulting classical evolution 
matrix needs to be multiplied by {v — 1)^^ to make it a legitimate (probability 
conserving) evolution operator. The resulting evolution does not permit back-scatter, 
and therefore, the secular equation fj43l) can be computed using Stark's C function ([8]), 
with Y = The product over the set of non back-scattering primitive periodic 

orbits becomes identical to the one appearing in the Ihara zeta function (jlj). Using (j5]), 
one finally gets, 

^M.(/i)=n(i-^,_i 

n ^ 

= fl-(^)^y"detf/(^)(l + ^)-^c) . 
\ V — 1 J \ V — 1 V — 1 J 

Thus, the spectrum mj of consists of r — 1 fold degenerate eigenvalues at m^-^"* = 
±;j^, and the rest which can be computed from the spectrum of the discrete Laplacian 
A. 



^ ^"^ (44) 



^± _ (^-A,)±V(.-A,)^-4(.-l) 
The eigenvalue of the Laplacian corresponds to the eigenvalues 1 and of MK The 
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gap in the classical evolution spectrum is determined by the first non zero eigenvalue of 
L. 

The comment above may have interesting and novel consequences going beyond its 
anecdotal appearance. In quantum graphs, one can choose vertex scattering matrices 
from a much larger variety than offered by the vertex scattering matrices f|T6|) . Thus, it 
is possible to construct vertex scattering matrices which do not scatter backwards, but 
with equal scattering probability to the other vertices. The unitarity is maintained by 
a proper choice of the phases of the scattering amplitudes [22]. In such cases, and for 
f-regular graphs, the classical analogues are identical with M" and ( l44l) is applicable. 
Working with such systems is particularly interesting because in quantum chaos, the 
gap between the eigenvalue 1 and the rest of the spectrum determines whether the 
spectrum of the U matrix (and hence of the Schrodinger operator) display the statistics 
predicted by Random Matrix Theory, in the limit of large graphs. For non back- 
scattering dynamics (jHj) reduces the problem to the study of the spectrum of the 
Laplacian. The behavior of the gap in the laplacian spectrum of large graphs is an 
important subject in the theory of discrete graphs and number theory, related amongst 
others to the Ramanujan conjecture [21]. A detailed discussion of this connection will 
take the present manuscript far afield, and it is deferred to a future publication. 
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